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ABSTRACT 

We study the 4-200 keV spectral and temporal behaviour of the low mass 
X-ray binary 4U 1820-30 with INTEGRAL during 2003-2005. This source as 
been observed in both the soft (banana) and hard (island) spectral states. A 
high energy tail, above 50 keV, in the hard state has been observed for the first 
time. This places the source in the category of X-ray bursters showing high- 
energy emission. The tail can be modeled as a soft power law component, with 
the photon index of ~2.4, on top of thermal Comptonization emission from a 
plasma with the electron temperature of kT e ~ 6 keV and optical depth of r ~ 4. 
Alternatively, but at a lower goodness of the fit, the hard-state broad band 
spectrum can be accounted for by emission from a hybrid, thermal- nonthermal, 
plasma. During this monitoring the source spent most of the time in the soft 
state, usual for this source, and the >4 keV spectra are represented by thermal 
Comptonization with kT e ~ 3 keV and r ~ 6-7. 

Subject headings: X-rays: binaries - X-rays: stars - stars: individual: 4U 1820-30 
- stars: neutron 



1. Introduction 

4U 1820-30 is a low mass X-ray binary seen at 0.66" from the centre of the globular clus- 
ter NGC 6624. It was the first identified source of type-I X-ray bursts (Grindlay et al. 1976). 
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Distance has been estimated as d ~ 5.8-8.0 kpc (e.g., Shaposhikov & Titarchuk 2004; Ku- 
ulkers et al. 2003). Kuulkers et al. (2003) found their determination of d = 7.6 ± 0.4 kpc 
compatible with the peak X-ray burst luminosity being a standard candle in their sam- 
ple of sources, and equal to the Eddington limit for He. The binary system consists of 
a He white dwarf with mass of O.O6-O.O8M (Rappaport et al. 1987) and a neutron star, 
with the mass estimated by Shaposhnikov & Titarchuk as ~ 1.3M Q , orbiting at short pe- 
riod of 11.4 minutes (Stella et al. 1987). In X-rays, 4U 1820-30 is classified as an atoll 
(Hasinger & van der Klis 1989). However, its flux variation between the soft (banana) to 
the hard (island) state are quasiperiodic at ~170 d (Priedhorsky & Terrell 1984; Simon 
2003; Wen et al. 2006), which quasiperiodicity has been proposed to be due to tidal effects 
of a more remote third star (Chou & Grindlay 2001; Zdziarski, Wen & Gierlihski 2006). 
X-ray bursts occur only at low flux levels (Cornelisse et al. 2003) and the kHz QPOs are 
correlated with the flux (Zhang et al. 1998), proving that the variability is due to an intrinsic 
change of the source luminosity (and the accretion rate) rather than being due to variable 
either obscuration or viewing angle of an anisotropic emitter. X-ray spectra of the source 
from RXTE and BeppoSAX were fitted by thermal Comptonization with or without an 
additional blackbody (Bloser et al. 2000; Sidoli et al. 2001; Migliari et al. 2004). 

In this paper, we report the first detection of X-ray emission above 50 keV from 4U 
1820-30 in a hard spectral state. We studied the source in the 4-200 keV energy range 
with the JEM-X (Lund et al. 2003) and IBIS (Ubertini et al. 2003) instruments on board 
the INTEGRAL satellite (Winkler et al. 2003a), collected during the Galactic Centre Deep 
Exposure (GCDE) program (Winkler et al. 2003b). We combined three methods of analysis: 
Temporal, using the light curves in different energy bands; Photometric, with the Hardness- 
Instensity Diagram; and Spectral, by spectral fitting. 

2. Observations and Data Analysis 

The coded mask imager IBIS on board INTEGRAL has a field of view (FOV) of 19° x 19° 
or 9° x 9° when partially or fully coded, respectively. The X-ray monitor JEM-X (4-35 keV) 
has the fully coded FOV of 4.8° x 4.8°. An INTEGRAL revolution lasts 3 days, and consists 
of adjacent pointings of ~2000 s, called Science Windows (SCW; Winkler et al. 1999). 

We have analysed all the available data during which the source was within the IBIS/ISGRI 
and JEM-X fully coded FOV, so that the flux evaluation is not affected by calibration un- 
certainties in the off- axis response. This yields 308 and 51 SCWs for IBIS and JEM-X 
respectively, during revolutions 50-363. Data were extracted with the Off-Line Scientific 
Analysis (OSA) (Goldwurm et al. 2003) v. 5.1 software released by the INTEGRAL Science 
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Data Centre (Courvoisier et al. 2003). A 2% systematic errors to both JEM-X and ISGRI 
data sets has been added 1 . Fluxes have been normalized to ISGRI value. The spectral 
analysis was performed with the XSPEC package v. 11.3. 

3. Light curves and hardness intensity Diagrams 

Fig. 1 shows the light curves of the monitoring period (2003 March 12-2005 October 
5) from soft to hard X-rays, in the 4-10, 10-20 keV bands from JEM-X and 20-30, 30-60, 
60-120 keV with IBIS/ISGRI, as well as the corresponding RXTE/ ASM 2 light curve. The 
~170-d quasiperiodic variability can be seen in the ASM data. The INTEGRAL light curves 
are marked with vertical lines to identify five epochs (A, B, C, D, E) during which we have 
performed detailed spectral analysis (see Table 1). Epoch A corresponds to the joint IBIS 
and JEM-X count rate maximum (~530 mCrab at 4-10 keV). Epoch C corresponds to the 
minimum of the 4-10 keV count rate (~100 mCrab), but to high count rates from 20 to 
60 keV, implying a spectral hardening across the 4-60 keV band. We also studied epochs 
before (B), and after that event (D, E). 

Fig. 2 shows hardness-intensity diagrams with data from JEM-X (a) and IBIS / JEM-X 
(b) (for simultaneous pointing and when possible because of different FOV). In Figs. 2 we 
can see the banana state, forming nearly horizontal bands across the A, B, D, E epochs. 
Thus, the source evolution is mostly in the flux, with the 4-20 keV spectral shape nearly 
constant. The upper part of the red points (epoch C) in Fig. 2a corresponds to the island 
state. In this case the 4-20 keV flux is almost constant while the hardness shows large 
changes. This indicates spectral pivoting somewhere in this energy band. Fig. 2b allows us 
to study spectral variability of the source at the energy band 10-30 keV, grater than that 
usually studied. In this diagram, we see similar behaviour in the horizontal banana state, 
while there is only one SC W for the island state, with a large hard colour. These findings 
are confirmed by spectral analysis as discussed below. 

4. Spectral analysis 

Data of the five epochs (see Table 1) have been fitted with a number of models and 
their combinations, namely, thermal and hybrid Comptonization, blackbody, disk blackbody 
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and a power law, all absorbed by the interstellar medium at 7V H = 3 x 10 21 cm -2 (Bloser et 
al. 2000). First, a thermal Comptonization model, CompTT (Titarchuk 1994), sufficiently 
model the spectra of the banana states (A, B, D, E). We obtain the electron temperature, 
kT c ~ 2.7-3.1 keV, and optical depth r ~ 6-7, as reported in Table 2. This is consistent with 
the colour-intensity diagram, where we found the spectral changes mostly in flux (most likely 
proportional to the accretion rate) and only slightly in their shape. The temperatures of the 
seed photons (Table 2) have been frozen at their best-fit values, given that those photons 
are entirely below the fitted band. The data, model, and residuals for the spectrum A 
(characterized by highest flux and unabsorbed bolometric luminosity of 7.7x 10 37 (rf/5.8 kpc) 2 
erg s -1 ) are shown in Fig. 3a, b, and all unabsorbed models, in Fig. 4. An addition of a 
blackbody component does not significantly improve the fits (although its presence is hinted 
by the negative residuals at the lowest energies, Fig. 3b), and thus the confidence ranges of 
the temperature, kT^, are not constrained. Still, we find kT^ ~ 1.5-2.4 keV at the best 
fits, in agreement with previous study of Bloser et al. (2000). 

The spectrum of the island state (C, red) shows an evident excess above 50 keV that 
cannot be fitted by thermal Comptonization only, as shown in Fig. 3e (the resulted xl is 
1.6 with 59 d.o.f.). An addition of a soft (r ~ 2.4) power law improves the fit significantly 
(xl = 1-0, see Figs. 3c, d), with the chance probability from the F-test (Bevington & 
Robinson 1992) of 6.3 x 10~ 7 . The parameters of the thermal plasma are also different from 
those for the banana state, with a higher kT e (~ 6 keV), and a lower r (^4), see Table 2. 
Fig. 4 shows the unabsorbed spectrum with its two components. By adding a disk blackbody 
component (Mitsuda et al. 1984) to the thermal Compton model, we have also obtained a 
good fit, with xl — 1-0. However, the required inner temperature, kT in ~ 5.4 keV, is clearly 
unphysical. Whereas adding a single blackbody gave a worse fit, as well as an unphysically 
high temperature, ~ 4 keV. 

We have investigated whether the tail can be an artefact of calibration inaccuracies. 
Analysing the data with previous software version of the OSA, the statistical significance of 
tail remains unchanged. 

The island-state data (C) have been fitted also with the highly-accurate Comptonization 
model COMPPS (Poutanen & Svensson 1996). Spectra agree very well with results from 
Monte Carlo Comptonization (Zdziarski et al. 2000, 2003). We assumed here spherical 
geometry of the Comptonizing plasma. Also, we have taken into account Compton reflection 
(Magdziarz & Zdziarski 1995) detected in other atolls (e.g., Zdziarski, Lubihski & Smith 
1999; Barret et al. 2000; Gierlihski & Done 2002), as well as a disk blackbody. We have 
obtained a very good fit, x^(d-o.f.)= 0.92(54) for a purely-thermal Comptonizing plasma 
with kT c ~ 55 keV and r ~ 0.6, and the solid angle subtended by reflector of f2/27r ~ 0.69 
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(assuming the source inclination of i — 60°). The bolometric luminosity of this model 
is 1.8 x 10 37 (<i/5.8 kpc) 2 erg s _1 . However, similar to the problem with the model using 
CompTT (see above), the inner disk temperature is rather high, kT in ~ 2.4 keV, as well as 
the temperature of the seed blackbody photons upscattered in the Comptonizing plasma is 
unphysically high, kT ~ 3.9 keV. 

In order to constrain the model to plausible values of the blackbody temperatures, 
we have considered a model with kT allowed to be < 2 keV. Also, we have constrained 
Q/2tt < 1. The resulting purely-thermal model provides a worse fit, x^(d.o.f.)= 1.09(54), 
with kT e ~ 24 keV, r ~ 1.9, kT[ n ~ 1.4 keV, similar to the parameters found in the atoll 
source 4U 1608-52 (Gierlihski & Done 2002). To improve the fit of this model, in particular 
to account for the high-energy residuals (appearing as a high-energy tail), we have added 
nonthermal electrons to the thermal plasma, as done to account for some spectra of both 
black-hole (Gierlihski et al. 1999; Wardzihski et al. 2002) and neutron-star (Farinelli et al. 
2005) binaries. We have obtained xl(d-o.f.) = 1.07(52) for a model with the nonthermal 
tail with the power-law index of p ~ 1 above the electron Lorentz factor of 1.6 (and the 
Maxwellian distribution below it), and other parameters similar to that of the previous 
model. We stress, however, that this model is significantly more complex than the one with 
CompTT and a power law presented above, and still yields significant residuals at energies 
> 80 keV. 



5. Discussion and Conclusions 

We have presented the first detection of hard emission above 50 keV in the island state 
of the atoll source 4U 1820-30 using IBIS. Up to now, 4U 1820-30 has not been included 
in the sample of the high-energy emitting bursters (Bazzano et al. 2006) and indeed neither 
BeppoSAX, CGRO/BATSE nor RXTE detected this source above 50 keV (Bloser et al. 2000 
and references therein). The previous lack of detection might have been due to either the 
source spending most of the time in the soft state, the weakness of the emission above 50 
keV in the island state, or the strong emission at > 50 keV appearing in that state only 
occasionally. 

In our island-state observation, the relative flux contribution for energies > 60 keV is 
~10%. The ratio of the 4-20 keV flux to the 20-120 keV ones is ~2, whereas it is in the range 
20-27 for the soft (banana) state. We point out that we detected emission above 50 keV 
when the 1-20 keV luminosity was ~ 1.4 x 10 37 (rf/5.8 kpc) 2 erg s -1 , i.e., close or somewhat 
below the critical value at which the X-ray binaries become hard emitters as proposed by 
Barret et al. (1996, 2000). 
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While the origin of the thermal Comptonization component in the soft state appears 
well understood as emission from hot electrons that upscatter soft photons coming from the 
accretion disk and the neutron star surface, the origin of the high-energy tail is less clear. 

The radio emission in the X-ray binaries is suggested to originate from jet-like out- 
flows which may also produce X-rays (Fender 2004). In particular, there are indication that 
both black-hole candidate and neutron stars sytems, when in hard state, produce steady jet 
(Migliari & Fender 2006). Also, in the hard state, positive correlation between radio and 
X-ray fluxes has been reported (Migliari & Fender 2006). 4U 1820-30 is one of the few radio 
detected atoll source. The only radio emission reported so far was during the soft state 
(Migliari et al. 2004), but it is probable that radio emission exist also in the hard state, as 
e.g. for 4U 1728-34 (Migliari et al. 2003). Furthermore, the power-law tail detected by us is 
relatively soft, with T ~ 2.4, and dominates over the Comptonization emission already for 
E < 4 keV, as shown in Fig. 4. Integration to its minimum energy, E min , yields the luminos- 
ity in this component oc E^^. If this component is due to nonthermal synchrotron emission, 
E m i n corresponds to the self-absorption turnover, typically in the IR or optical range, as it is 
the case in another ultracompact binary, 4U 0614+091 (Migliari et al. 2006). Then the tail 
luminosity would be ~10 times that of X-rays (see also discussion in Zdziarski et al. 2003), 
making the island state more luminous than the banana state, with the maximum of the 
overall EFe at ~1 eV, neither of those seen in X-ray binaries. 

Also, the presence of the jets or outflow could be due to magnetic reconnection avents 
above the accretion disk (Shibata 2005). 

A more likely origin of the power law tail appears to be emission of nonthermal elec- 
trons in a hybrid (i.e., containing also thermal electrons) plasma, as discussed in Section 
4. Nonthermal tails have been commonly seen in soft states of black- hole binaries (e.g., 
Cyg X-l, Gierlihski et al. 1999) and neutron-star Z-sources, e.g., GX 17+2 (Di Salvo et al. 
2000; Farinelli et al. 2005), Sco X-l (D'Amico et al. 2001), GX 349+2 (Di Salvo et al. 2001), 
GX 5-1 (Asai et al. 1994), while they are difficult to detect in the atoll sources, e.g., 4U 
0614+091 (Piraino et al. 1999). On the other hand, nonthermal tails have been found in the 
hard state of some black-hole binaries, e.g. Cyg X-l (McConnell et al. 2002) and GX 339-4 
(Wardzihski et al. 2002). 

The 4U 1820-30 high-energy part of the spectrum (E > 10 keV) can be well fitted with 
purely thermal Comptonization at kT c ~ 50 keV, as discussed in Section 4. This electron 
temperature is higher than that seen in 4U 1608-52 (Gierlihski & Done 2002), but the 
Comptonization temperature of atolls in the island state has not yet been determined for a 
sample of sources, and kT e ~ 50 keV is similar to that seen in the hard state of black-hole 
binaries (e.g., Zdziarski & Gierlihski 2004). Then, the problem with fitting the spectrum 



-7- 



at < 10 keV with this model may be due to the complexity of the X-ray spectra in the 
hard state, including contributions from the disk, the corona above the disk, optically-thin 
and optically-thick parts of the boundary layer and the uncovered part of the neutron star 
surface. 

A different thermal Comptonization model has been considered by Markoff et al. (2005). 
They proposed thermal synchrotron self-Compton emission in a jet as a model for the hard 
state of black hole binaries. However, their model requires an extremely high electron tem- 
perature, kT e ~ 3-4 MeV. This would result in the thermal spectral cutoff in black hole 
binaries at energies > 3kT e ~ 10 MeV, which has not been seen in any object yet. Although 
we cannot formally rule out a contribution from this process, postulating the presence of ~3 
MeV plasma in 4U 1820-30 appears rather speculative at this time. 

In summary, we have discovered strong high-energy emission from 4U 1820-30. The 
high energy emission is harder than that from some other atolls in the island state. Its most 
likely interpretation is either emission of some nonthermal electrons or thermal emission from 
plasmas with a relatively high temperature, kT e > 50 keV. 

This research has made use of data obtained through the INTEGRAL Science Data 
Centre, Versoix, Switzerland. The authors thanks M. Federici for the continuous effort to 
update the INTEGRAL archive and software in Rome and G. De Cesare, L. Natalucci and 
F. Capitanio for the scientific and data analysis support. Particular thanks are due to P. 
Goldoni for making his private data available for analysis. This work has been supported by 
the Italian Space Agency through the grant I/R/046/04. AAZ has been supported by the 
Polish grants 1P03D01827, 1P03D01128 and 4T12E04727. 
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Table 1: The log of the data used for spectral fitting. The indicated colours are used in the 
light curves, Fig. 1, and hardness-intensity plots, Fig. 2. 



data set 


Time (MJD) 
start— stop 


rev. # 
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IBIS ; JEM-X a 


Obs. type 


A (blue) 
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GCDE 


B (green) 
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GCDE 


C (red) 
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GCDE 



Q JEM-X2 for the A obs. and JEM-X1 for the B-E obs. 
b A03 target-of-opportunity observation, PI P. Goldoni. 
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Fig. 1. — The 2003-2005 light curves of 4U 1820-30, with the detector count rate in each 
band given with respect to the corresponding Crab count rate. The panels are marked with 
the energy range and detector, and present the INTEGRAL SCW data points except for the 
top panel, which gives the RXTE/ ASM 1-day averages. The lines mark the data sets (A, B, 
C, D, E) used for joint IBIS and JEM-X spectral fits (see Table 1). [See on-line edition for 
the coloured version of this figure]. 
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Fig. 2. — (a) Hard colour-intensity diagram; each point corresponds to 100 s. (b) Very hard 
colour-intensity diagram; each point corresponds to a SCW. The colours identify the A-E 
data sets, see Table 1 and Fig. 1. The typical error bars are shown at the bottom right 
corners. 



Table 2: Spectral fitting results for the JEM-X and IBIS broad-band spectra. The model is 
CompTT for A, B, D, E, and a CompTT + a power law for C. 



parameters 


A (blue) 


B (green) 


C (red) 


D (purple) 


E (orange) 


kT (keV) a 


0.2 


0.4 


1.5 


0.3 


0.2 


kT e (keV) 
r 

nornicompTT 


o 71 +0.08 
z - ' 1 -0.07 

7 04 +a32 

' - u ^-0.29 
9 o 4 +0.21 
z -°^-0.19 


q 19+O.IO 

°- ±z -o.io 

K QC+0.24 
J.OO_o 23 

u - do -0.06 


6 10 +0 - 77 

oq+0.79 
°- oy -0.65 

3.79j£$ xlO" 2 


2 qn +0 - 20 

z - yu -0.18 

1 06+ - 21 

i - uu -0.18 


9 OO+0.07 

z - oo -0.06 

6 71 +0 - 24 

°- ' i -0.22 
1 1-7+O.I2 
L -°' -0.11 


r 

norm p ; 






9 qcr +0.10 
-0.11 

60+ - 22 

U.UU_ 24 






Xl(d.0l) 


1.13(44) 


0.96(46) 


1.0(57) 


1.07(37) 


1.14(43) 


-F4-20keV (erg s" 1 cm" 2 ) 


9.1xl0" 9 


3.8xl0" 9 


2.1xl0" 9 


4.2xl0" 9 


5.8xl0" 9 


i^o-eokev (erg s" 1 cm" 2 ) 


3.3xl0" 10 


1.8xl0" 10 


7.4xl0" 10 


2.2xl0" 10 


2.5xl0" 10 


^60-i20keV (erg s" 1 cm" 2 ) 






1.5xl0" 10 







"Fixed parameters 
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Fig. 3. — Data, models and residuals of the soft spectrum A (a, b) and the hard spectrum 
C (c, d) fitted by thermal Comptonization without and with, respectively, a power law 
component. The panel e shows the residuals for the spectrum C fitted without the power 
law component. The high energy excess at > 50 keV is clearly visible. 
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Fig. 4. — The unabsorbed best-fit models of the A-E data sets, see Table 1. The model for 
the green, purple, orange and blue spectra is thermal Comptonization, while a power law 
(dotted curve) is added in the red spectrum (C). Note that the power law dominates at low 
energies, and thus it is probably only a phenomelogical, but unphysical, description of the 
spectrum. 



